Coupling between magnetic ordering and structural instabilities in perovskite 

biferroics: A first-principles study 
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We use first-principles density functional theory-based calculations to investigate structural insta- 
bilities in the high symmetry cubic perovskite structure of rare-earth (R = La, Y, Lu) and Bi-based 
biferroic chromites, focusing on V and R point phonons of states with para-, ferro-, and antiferro- 
magnetic ordering. We find that (a) the structure with G-type antiferromagnetic ordering is most 
stable, (b) the most dominant structural instabilities in these oxides are the ones associated with ro- 
tations of oxygen octahedra, and (c) structural instabilities involving changes in Cr-O-Cr bond angle 
depend sensitively on the changes in magnetic ordering. The dependence of structural instabilities 
on magnetic ordering can be understood in terms of how super-exchange interactions depend on the 
Cr-O-Cr bond angles and Cr-O bond lengths. We demonstrate how adequate buckling of Cr-O-Cr 
chains can favour ferromagnetism. Born effective charges (BEC) calculated using the Berry phase 
expression are found to be anomalously large for the A-cations, indicating their chemical relevance 
to ferroelectric distortions. 



I. INTRODUCTION 

A ferroic is a material which exhibits spontaneous and 
switchable ordering of electric polarization or magnetiza- 
tion or elastic strain. Materials exhibiting more than one 
of such orderings termed 'multiferroics'- have recently be- 
come the focus of much research^. Most of the biferroics 
investigated in recent years are ABO3 oxides with per- 
ovskite structure. The d°— ness or the zero occupancy 
of tranistion metal B cation is known chemically to fa- 
vor ferroelectricity 2 -. Hence, the availability of transition 
metal <i-electrons in the perovskite oxides necessary for 
magnetism, reduces the tendency for off-centering ferro- 
electric distortions 2 making multiferroics relatively rare. 
How the ordering of d— electronic spins of the B cation in- 
fluence ferroelectric or other competing structural insta- 
bilities has not yet been explored and understood. The 
coupling between magnetic ordering and structural in- 
stabilities is expected to involve interesting physics and 
is of direct relevance to technological applications 3 such 
as multiple state memory elements and novel memory 
media. 

There are at least three families of high tempera- 
ture biferroic materials. Bi- based perovskite oxides, 
like BiMnO-A BiCrOa- 5 , and BiFeO^i are known to 
be promising biferroics. Ferroelectricity in these materi- 
als arises from the stereochemical activity of the 6s lone 
pair electrons of Bi. Hexagonal rare earth manganates 
LnMnOa and InMnOs are biferroics which exhibit im- 
proper or geometric ferroelectricit y 8 ' 9 ' 10 ' 11 . Rare earth 
chromites LnCrO ^ 12 ' 13 (with Ln = Ho, Er, Tm, Yb, Lu 
or Y) have been recently shown to be biferroic; YCr03 
has been showni 2 - to exhibit canted antiferromagnetic be- 
havior below 140 K and a ferroelectric transition around 
473 K. Similarly, LuCrOs becomes a canted antifcrro- 
magnet below 115 K, and is ferroelectric below 488 K-i 3 -. 
The absence of any ferroelctricity in LaCrC>3, has been 



attributed to the large size of the La 3+ ion, in comparison 
withY 3 +. 

However, small values of polarization reported for 
these materials (about 2/iC/cm 2 for YCr03^2 and 
6/iC/cm 2 of BiFeOA) inspite of large A-cation off- 
centering distortions remain a puzzle. More recently", 
a new concept of 'local non-centrosymmetry' in YCrC>3 
has been proposed to account for the small value of 
polarization observed. Perovskite oxides are known to 
have many competing structural instabilities^ and this 
competition is further enriched by the magnetic insta- 
bilities. The coupling and competition between various 
magnetic and structural ordering can be partly respon- 
sible for weak ferroelectricity or the possibility of local 
non-centrosymmetry. Our goal here is to investigate this 
issue through determination of various structural insta- 
bilities for different magnetic orderings, with a focus on 
rare earth chromites. 

We present results of detailed electronic structure and 
frozen phonon (at the Y and R points) calculations for a 
set of five materials (LuCrC>3, YCr03, LaCr03, BiCrOs 
and YFeOs) in the cubic phase with three different mag- 
netic orderings (para-, ferro- and antiferromagnetic). In 
Section II, we briefly describe the methods used in cal- 
culations here. In section III, we report results for struc- 
tural energetics, the electronic density of states (DOS) 
and the Born Effective charges (BEC) for the cubic phase 
of these materials with different magnetic orderings. In 
Section IV, we report results for structural instabilities 
in chromites and compare them with those in a related 
compound YFe03. Since the many-electron correlations 
are important in magnetic oxides, we estimate their ef- 
fects on the structural instabilities through use of the 
Hubbard parameter U 16 . Using our results for structural 
instabilities, we show how certain nonpolar structural 
instabilities can cooperatively stabilize ferromagnetism. 
Our work reveals how the structural instabilities of these 
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biferroic oxides depend on the size of A-cation (with the 
same B-cation), and a change in B-cation. We interpret 
the results using arguments based on superexchange 17 , 
and the well-known Goodenough-Kanamori rulesi^. Fi- 
nally, we summarize in Section V. 
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FIG. 1: Perovskite structure: Cell doubled along the (111) 
direction, to represent G-type antiferroamgnetic ordering. 



II. METHODOLOGY 

Our calculations are based on first-principles pseu- 
dopotential based Density Functional Theory within a 
generalized gradient approximation (GGA)i^ as imple- 
mented in the PWSCF package 2 ^. The interaction be- 
tween ions and electrons is approximated with ultra- 
soft pseudopotentials^i treating explicitly 11 electrons 
[(n-l)s 2 (n-l)p 6 (n-l)d 1 and ns 2 ] in the valence shell of 
Lu(n=6), La(n=6) and Y(n=5). We consider 6 valence 
electrons for Oxygen [2s 2 2p 4 ] and 14 electrons for Cr 
[3s 2 3p 6 3d 5 4s 1 ]. We used a plane wave basis with kinetic 
energy cut off of 25 Ryd (150 Ryd) to represent wave 
functions (density). For cubic structures, we sample the 
Brillouin zone using a 5 x 5 x 5 Monkhorst Pack Mesh^ 2 -, 
and a denser mesh (6x6x6) and higher energy cut-off 
(30 Ryd) for energy differences. Phonon frequencies cal- 
culated with these larger parameters do not differ much 
from those calculated with a lower cut-off. We perform 
spin polarised calculations by initializing different spins 
on neighbouring magnetic ions; For paramagnetic order- 
ing we initialise a zero spin on Cr ions. To represent anti- 
ferromagnetic ordering, the unit cell is doubled along the 
(111) direction (see Fi^TJ). To investigate structural in- 
stabilities in the prototype cubic structure, we determine 
its dynamical matrix using frozen phonon calculations. 
We use a finite difference form of the first derivative to 
compute an element of the force constant matrix: 



TABLE I: Lattice constants of various oxides in the cubic 
perovskite structure with experimental unit cell volumes. 





Stress(GPa) 


a (A) 


PM 


FM 


AFM 


LuCr0 3 


3.77 


-4.712 


3.302 


1.429 


YCr0 3 


3.79 


-3.875 


3.883 


2.240 


LaCrOa 


3.88 


-7.043 


-0.083 


-1.564 


YFeOa 


3.83° 


-15.8 


-5.9 


2.5 


BiCr0 3 


3.85 6 


-15.4 


-9.7 


-11.2 
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where Fi a is the Hellman-Feynman force acting on the 
t h atom in a direction, and, Ujg the displacement of the 
/ atom in (3 direction with respect to the equilibrium 
structure. We used A=0.04 A, about 1% of the lattice 
constant. The dynamical matrix is then calculated from 
the force constant matrix, 
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whose eigenvalues correspond to the square of the phonon 
frequencies (w 2 ). 

In periodic systems, the dynamical charge tensor or 
Born effective charge tensor can be defined 23 as the co- 
efficient of proportionality between the macroscopic po- 
larization created in direction (3 and a rigid displacement 
of the sublattice of atoms j in direction a, 
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fl being the unit cell volume. The polarization is deter- 
mined using the berry phase formalism 2 ^ as implemented 
in the PWSCF package. 



III. PROPERTIES OF THE CUBIC 
PEROVSKITE STRUCTURE 

We have determined electronic structure of RCr03 
compounds (R=Y, Lu, La) in the high symmetry cu- 
bic structure with different magnetic orderings. This is 
accomplished through calculations with different initial 
guesses for atomic spin polarization and optimizing with 
respect to spin density. All our calculations are for the ex- 
perimental unit cell volumes, as ferroelectricity is known 
to be sensitive to lattice constants or pressure (see lat- 
tice constants listed in Table [I]). In many magnetic com- 
pounds, a change of magnetic ordering causes a stress 
which induces a structural distortion 2 ^. This concept of 
'magnetic stress' was introduced to describe structural 
phase transitions that are induced by magnetic order- 
ing, and applied to materials with degenerate (usual e 9 ) 
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orbitals 2 ^. In these materials as well, a change in mag- 
netic ordering (with fixed lattice parameters) produces a 
change in stress (see Table QJ. For LuCrC>3 and YCrC>3, 
the introduction of spin polarization produces a change 
from compressive stress in the paramagnetic phase to a 
tensile stress, with stress being minimum in the antifer- 
romagnetic structure (the lowest energy ordering). In 
LaCrOs, the stress remains compressive with all three 
magnetic orderings. 



Electronic structure of YCrOa, 
LaCrOa 



LuCrOs and 



1. Cubic Paramagnetic(PM) Structure 
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FIG. 2: Density of States for cubic PM YCr0 3 , LuCrOa and 
LaCrOs. The fermi level (indicated by a dashed line) has 
been set to zero in all the three cases. 

First, we present results for the highest symmetry cu- 
bic structure with no spin polarization. Although this 
state is experimentally inaccessible, it provides a useful 
reference for understanding the spin-polarized structures 
discussed later in the paper. The plotted energy range 
is from -10 to 4 eV, and the lower lying semicore states 
have been omitted for clarity. In the PM cubic YCrC>3, 
LuCr03 and LaCrC>3, (see Fig [5]) there is high density 
of electronic states at the Fermi level, driving the system 
towards a Stoner instability 28 . This suggests that this 
phase should be unstable with respect to spin polariza- 
tion and/or structural distortions. The contribution of 
various orbitals to the DOS can be understood better by 
examining the orbital projected density of states (see Fig 
[3]) which show that, the contribution between -8 to -3 eV 
is mainly from Oxygen 2p orbitals. Cr d-orbitals con- 
tribute predominantly to the peaks at the Fermi level. 
In contrast, the contribution from the Lu d-orbitals is 
substantial only 2eV above the Fermi level. 



2. Cubic Ferromagnetic(FM) Structure 

Ferromagnetic cubic structure is simulated by initial- 
izing spins on both Cr ions in the same direction. In 



u 1 

i 

c 



I 

I 

LuW a 

Lubl 1^ 



Da, 

' >\ r" 

j i y r 




Ci3flw.. 





-1ft -ft 



Energy <evj 



FIG. 3: Orbital resolved density of states for cubic PM 
LuCrOs. The high density of states at the Fermi level hints 
that it is an unstable phase. 
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FIG. 4: Total and orbital resolved Density of states for cubic 
ferromagnetic YCr03, LuCrOs and LaCrOs with the Fermi 
level set to zero in all the three cases. 



all chromites studied here, the ferromagnetic structures 
have a magnetic moment of 3 /is in accordance with the 
Hund's rule value expected for a d 3 configuration. The 
majority spins are represented by the solid line on the 
positive Y axis, and the minority spins on the negative 
Y axis. The introduction of spin polarization reduces 
the energy by approximately 2 eV per unit cell. The 
source of stabilization is clear from the density of states 
(see Fig 0} which reveals opening of a gap at the Fermi 
level. The states corresponding to non-magnetic atoms 
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are unchanged in comparison with PM ordering. The 
down-spin Cr 3d states are split off from the O 2p states 
creating a wide gap for the minority states. The up-spin 
Cr 3d states hybridize with the O 2p states and there 
is a very small gap for the majority carriers. The den- 
sity of electronic states at the Fermi level is still finite 
having a small contribution from the Cr d-orbitals. This 
hints that either the ferromagnetic phase may not be the 
most stable, and that either an antiferromagnetic (AFM) 
spin arrangement could lower the energy of the system, 
or that the cubic structure is unstable and a structural 
distortion will lower the energy of the system. Since Cr 3+ 
is a d 3 ion, it is Jahn- Teller inactive, and the structural 
distortions (if any) probably involve the A-cation (at the 
corners), or the oxygen anions. 



respect to PM ordering is much more (around 2eV) than 
the gain in energy in going from the FM to the AFM 
structure (see Table UTT|) . LuCrC>3, like YCr03, is also 
found to be insulating with the introduction of a gap at 
the Fermi level (see Fig. [5|. Both spin channels have 
identical density of states consistent with an AFM spin 
arrangement. 

From the orbital projected density of states for AFM 
arrangement, we find that the t2 9 orbitals of Cr are fully 
occupied (and constitute the HOMO), whereas the e g 
orbitals are unoccupied. Although the LUMO consists 
of Cr d-orbitals (about 2eV above the Fermi level), the 
Lu d states also appear within the same energy range 
(see Fig. [5]). 



3. Antiferromagnetic Structure 

We simulated antiferromagnetic structure by initializ- 
ing antiparallel spins on the two Cr ions in the super- 
cell. It is well known that, the superexchange between 
e g orbitals of adjacent ions connected through oxygen 
with a 180° metal-oxygen-metal bond angle, is much 
stronger than the interaction between the correspond- 
ing t2 g orbitals, since the former is mediated by stronger 
dpcr bonds as compared to the weaker dp7r bonds in the 
latte r 17 i 18 . So, we expect a superexchange interaction in 
which there is a weak coupling between the t2 S orbitals 
of the adjacent Cr atoms giving rise to an antiferromag- 
netic interaction. Further, this coupling will be stronger 
in the cubic structure as the bond angle between Cr-O-Cr 
is 180°, as compared to a distorted structure. From su- 
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FIG. 5: Density of states for antiferromagnetic YCrC>3, 
LuCrO-3 and LaCr03. 

perexchange arguments applied to d 3 configurations, the 
structure with G-type antiferromagnetic ordering hav- 
ing rhombohedral symmetry should be most stable. We 
consider collinear spins assuming that the canting of the 
spins would be small. We find that the AFM structure 
is lower in energy by about 0.4 eV than the FM phase. 
We note that this gain in energy by FM ordering with 
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FIG. 6: Orbital resolved density of states for antiferromag- 
netic LuCrOs. 



B. Born Effective Charge (BEC) 

The effective charge tensors have been calculated from 
polarization differences between the perfect and distorted 
structures in the AFM phase. The anomalous values of 
Z* so obtained indicate that a large force is felt by a given 
ion due to small macroscopic electric field, thus favoring a 
tendency for off-centering and toward a polarized ground 
state. The effective charge of the A-cation (see Table 
HI)) is about the same for the three cases. For LaCrC>3, 
we find a larger BEC on Cr and one of the oxygen atom 
moving along the bond. This is possibly because of larger 
Cr-0 bond length (arising from larger size of La cation, 
see Table IIIII) and correspondingly greater contribution 
from the long-range charge transfer. We expect from this, 
that the superexchange interaction in LaCrOa should be 
stronger as well. 



5 



TABLE II: The XX component of Born effective charge tensor 
for AFM LaCrC>3, LuCrOa and YCr0 3 . Nominal charges are 
indicated in brackets. 













LuCr0 3 


4.42(3) 


3.43(3) 


-2.56(-2) 


-2.62(-2)) 


YCrOa 


4.45(3) 


3.44(3) 


-2.62(-2) 


-2.66(-2)) 


LaCr0 3 


4.5(3) 


3.76(3) 


-3.82(-2) 


-2.22(-2)) 



TABLE III: Relative Energies of different Magnetic phases, 
Cr-0 bond lengths, and Neel's temperatures for cubic 
LuCr0 3 , YCr0 3 and LaCr0 3 (Energy of the PM phase has 
been set to zero). 





PM 


FM 


G-AFM 


Cr-O Bond length 


Tat 


LuCrO a 


0.0 


1.89 eV 


2.3 eV 


1.88(A) 


115 K a 


YCrOa 


0.0 


2.02 eV 


2.4 eV 


1.90(A) 


140 K° 


LaCr0 3 


0.0 


2.04 eV 


2.5 eV 


1.94(A) 


282 K b 



a Neel's temperature taken from Ref. 1 1-JI I il . 
^Ref.Q. 



IV. STRUCTURAL INSTABILITIES 
A. Coupling with Magnetic Ordering 

In order to represent G-type antiferromagnetic order- 
ing which has been shown to be most favourable energet- 
ically, we use a unit cell doubled along the (111) direc- 
tion. We determine structural instabilities in this struc- 
ture with different magnetic orderings. A single unit cell 
has 10 atoms which results in 30 phonon branches: 3 ac- 
coustic (which have zero frequency at k=(0,0,0)) and 27 
optical, some of which are triply degenerate. We are in- 
terested mainly in optical modes with imaginary phonon 
frequencies corresponding to instabilities in the structure. 
Doubling the unit cell along the (111) direction, gives us 
access to zone boundary phonon modes (R-point) which 
form the dominant structural instabilities in this struc- 
ture, along with the zone-center modes. 

In the paramagnetic phase, both YCrOs and LuCrOs 
exhibit a zone center instability at 116.5 and 144.8 cm -1 
respectively, which is a polar mode (with Tis symme- 
try) involving mainly the off centering of A-cation. This 
instability, however, is absent in LaCrOs. We find two 
instabilities in the FM phase: and 1^5 modes. The 
non-polar mode involves oxygen displacements only, 
and is strongly unstable in the FM phase. The Tis mode 
involves the A-cation (rare earth ion) moving in a direc- 
tion opposite to that of the oxygen cage and Cr atom 
resulting in a ferroelectric polar structural distortion. 
Note that the Cr atom moves in the same direction as 
the oxygen ion, in contrast to the behavior of Ti ion in 
BaTiO^ and PbTiOiA but similar to the behavior of 
Mn in BiMnOgii and Cr in BiCrOgl 

With AFM spin arrangement, for LuCrC>3, we find 
three triply degenerate instabilities, at 339.5 cm -1 , 145 



cm' 1 and the weakest at around 60 cm -1 . For YCr03, 
the corresponding instabilities are at 309 cm -1 , 140 cm' 1 
and 80 cm' 1 respectively. The strongest instability (at 
around 300 cm' 1 for the two materials) has R25 symme- 
try and corresponds to rotation of the corner connected 
oxygen octahedra. The next instability is the ferroelec- 
tric Ti5 mode (around 140 cm' 1 ). The weakest insta- 
bility (60 cm' 1 for LuCr0 3 and 80 cm' 1 for YCr0 3 ) 
has R15 symmetry and involves displacement of the A- 
cations (Lu and Y for our case) and small oxygen dis- 
placements; these are antiparallel in neighbouring unit 
cells. In LaCrOs, we find only two triply degenerate in- 
stabilities. The first instability at around 220 cm' 1 cor- 
responding to the oxygen rotations (the R25 mode) and 
the second close to 18 cm' 1 having Tis symmetry (see 

Fig ED. 
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FIG. 7: Eigenvectors of the unstable P point phonon modes: 
Ti5 and T25 modes 

For the rare earth chromites studied in this paper, un- 
stable R-point modes in the AFM phase do not change 
significantly with magnetic ordering. T point instabili- 
ties in contrast depend strongly on magnetic ordering. 
Only the high frequency phonons are affected in going 
from para- to antiferromagnetic phase and all T15 and 
i?i5 modes are softer in the FM phase. 

The T25 mode which brings about a significant change 
in the Cr-O-Cr as well as O-Cr-0 bond angle shows a 
spectacular change with magnetic ordering. This mode, 
highly unstable in the ferromagnetic phase (at around 
200 cm _1 ) becomes stable at around 50 cm" 1 with 
para- and antiferromagnetic orderings. This behavior, 
although not as pronounced, is also seen for this mode 
in LaCrOs. The R25 instability, also involving a change 
in the Cr-O-Cr bond angle, is not affected by the change 
in magnetic ordering possibly because the O-Cr-0 bond 
angle still remains unchanged. 

Another significant change is observed for the R'25 
mode, which involves a movement of the two B-cations 
(Cr in our case) in opposite directions. After introduction 
of spin-polarization, the R25' mode (close to 250 cm -1 for 
PM phase) becomes more stable at around 400 cm -1 for 
the FM and AFM phases. 

We compare our results with the Bi-based biferroic 
chromite, BiCrOs and find a similar behavior of the T25 
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mode here as well. We thus attribute this behavior to 
the B cation (Cr ion for the chromites) and expect it to 
be their general behavior. To interpret the general trend 
in phonon frequencies (see Fig. [5]), the following rules 
apply: 

1. The modes which involve a change in the Cr-O- 
Cr bond angle (as well as O-Cr-0 bond angle) are more 
stable in the AFM phases and are relatively less, or un- 
stable, in the FM phase. This behaviour is seen in the 
r 2 5 and r 15 modes for LuCrC^. 

2. Secondly, the modes which involve a change in the 
Cr-0 bond length tend to harden with the introduction 
of spin polarization, as observed for the R'25 and R'12 
modes. 
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FIG. 9: Changes in phonon frequencies with changes in lo- 
cation: YCr0 3 and YFe0 3 . 

crossover, at a displacement of 0.407 A, beyond which the 
FM phase is energetically favoured. The Cr-O-Cr bond 
angle, at the crossover point is found to be 153°, which is 
significantly different from, the value suggested by Good- 
enough and Kanamori for ferromagnetic superexchange 
interaction (130°) in 1951. On examining the Density of 
states beyond the crossover point (see Fig. [TTj) . we find 
more significant hybridization between Cr and oxygens 
for the FM phase. Secondly, the FM phase so stabilized 
is found to be insulating. 



FIG. 8: Changes in phonon frequencies of cubic LuCrOa with 
different magnetic orderings. Insets show similar curves for 
LaCrOa and BiCrOa. 

In order to study the effect of change in B-cation (see 
Fig©, we compare instabilities in YFeC>3 with YCrC>3. 
We find that a change in magnetic ordering has an op- 
posite effect on the unstable modes involving a change 
in Fe-O-Fe bond angle. A spectacular change is the R'2 
mode, an oxygen breathing mode, which softens in the 
FM phase as compared to the PM and AFM phases. The 
mode also shows a different behaviour, showing a sta- 
bilization with FM ordering. These differences are due 
to the filled e g orbitals in Fe 3+ , which are unoccupied in 
Cr 3+ , leading to a much stronger superexchange inter- 
action mediated by the e g orbitals, and have a different 
geometry dependence. 



B. Stabilization of Ferromagnetic Ordering 
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FIG. 10: Energy vs displacement corresponding to the 1^5 
mode for FM and AFM orderings in LuCrOs. FM becomes 
more stable than the AFM state for rhombohedral distortions 
greater than 0.407 A. 



The spectacular change in the frequency of the T25 
mode with FM ordering, prompts us to discuss whether 
ferromagnetism can be stabilized in these chromites by 
varying the Cr-O-Cr bond angle. Since the T25 mode is 
unstable only in the FM phase, we want to study the 
effect of freezing in a distortion of this mode. We have 
calculated the total energy as a function of T 2 5 displace- 
ments in the (111) and (100) directions, for the FM and 
AFM phases (see Fig [TU)) . For LuCrOa, we observe a 



C. Effect of Correlations 

As mentioned earlier in the paper, the LDA+U method 
has been successfully applied to describe the electronic 
structure of sytems containing localized d and / elec- 
trons where LDA sometimes leads to incorrect results^, 
and recently it has been applied to obtain structural pa- 
rameters that are in better agreement with experimental 
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FIG. 11: Density of states for FM and AFM phases at r 25 
displacement of 0.407 A (see Fig. \TU\ where the FM state is 
more stable). 



results than LDA or GGAii2£. In this work, we use a 
value of U = 3.0 eV, adapted from work on full struc- 
tural optimization of BiCrOgi. With the introduction 
of correlations through U parameter, we find that the 
modes in FM phase do not change significantly. Our re- 
sults for YCrC>3 (see Fig [12]) bear that only r 2 5, R'2 and 
R'12 modes are noticeably affected. Correlations lead to 
softening of the R'2 and R'12 modes by 50 cm -1 , and 
tend to harden the T25 mode in the PM phase. In an 
AFM spin arrangement, these R-point modes are hard- 
ened by approximately 20 cm -1 and the T25 instability 
is softened. 
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FIG. 12: Effect of correlations on the phonon modes of YCrOs 
with different magnetic orderings. For each mode, data on the 
left of the vertical dashed line represents estimates with GGA 
and the data on the right represents estimates with GGA+U 
(U=3.0 eV). 



V. SUMMARY 

In conclusion, we have determined structural instabil- 
ities of LaCrOs, LuCr03, YCr03 and BiCrOs in their 
cubic perovskite structures with different magnetic or- 
derings. Our finding that the G-type antiferromagnetic 
ordering is most stable can be explained with superex- 
change arguments. Ferroelectric structural instabilities 
in the cubic structures involve A-cation (Lu or Y) dis- 
placements, as indicated by the eigenvectors of the fer- 
roelectric Ti5 modes and an anomalous BEC of the A- 
cation. We find that certain phonon frequencies de- 
pend sensitively on magnetic ordering: the modes in- 
volving a change in bond-angle are stable (harder) with 
the antiferro- and paramagnetic ordering than in the FM 
state; on the other hand, the modes involving a change in 
Cr-0 bond length are softer in the paramagnetic phase 
and comparable in the FM and AFM states. The T25 
oxygen mode brings about a significant change in the Cr- 
O-Cr bond angle and is highly unstable in the FM phase, 
and corresponding structural distortion leads to stabiliza- 
tion of ferromagnetic ordering in these chromites. Among 
the competing structural instabilities the antiferrodis- 
tortive instability (R25 mode) is the strongest. Electron 
correlations are found to have little effect on the unstable 
phonon modes, but result in a slight change in a few of 
the stable phonon modes in the PM phase. We note that 
the effects of magnetic ordering on structural instabilities 
are quite different (in fact, opposite sometimes) in YFe03 
with respect to chromites, and a more detailed study is 
required to understand such couplings in ferrites. Origin 
of small polarization and/or local non-centrosymmctry 14 
is probably from the relatively weak ferroelectric instabil- 
ities and their competition with various structural mag- 
netic instabilities, and our work should be useful in for- 
mulating a phenomenological analysis of the same. 
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